Existence of Fermion Zero Modes and Deconfinement of Spinons in Quantum 
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We investigate the quantum antiferromagnetism arising from algebraic spin liquid via spontaneous 
chiral symmetry breaking. We claim that in the antiferromagnet massive Dirac spinons can appear to 
make broad continuum spectrum at high energies in inelastic neutron scattering. The mechanism of 
spinon deconfinement results from the existence of fermion zero modes in single monopole potentials. 
Neel vectors can make a skyrmion configuration around a magnetic monopole of compact U(l) gauge 
fields. Remarkably, in the monopole-skyrmion composite potential the Dirac fermion is shown to 
have a zero mode. The emergence of the fermion zero mode forbids the condensation of monopoles, 
resulting in deconfinement of Dirac spinons in the quantum antiferromagnet. 
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I. INTRODUCTION 

High T c superconductivity is believed to result from 
hole doping to an antifcrromagnctic Mott insulator. Hole 
doping to an antiferromagnetic Mott insulator destroys 
antiferromagnetic long range order, resulting in one 
quantum disordered paramagnetic Mott insulator that is 
considered to be the pseudogap phase in high T c cuprates. 
High T c superconductivity is expected to arise from fur- 
ther hole doping to the paramagnetic Mott insulator 
In this respect it should be a starting point for the theory 
of high T c superconductivity to understand the nature of 
the antiferromagnetic Mott insulator. 

Recently, the paramagnetic Mott insulator was pro- 
posed to be algebraic spin liquid, where spin 1 antifer- 
romagnetic fluctuations break up into more elementary 
spin 1/2 fractionalized excitations called spm<ms[lj. If 
the algebraic spin liquid correctly describes spin degrees 
of freedom in the pseudogap phase, its parent antifer- 
romagnet may include the trace of fractionalized spinon 
excitations. In the present paper we find a trace of de- 
confincd spinons in the quantum antiferromagnet. 

There are some experimental reports that the spin 
spectrum in the antiferromagnetic Mott insulator is dif- 
ficult to understand only by antiferromagnons 2J. Al- 
though the dispersing peaks observed in inelastic neutron 
scattering measurements can be interpreted as antiferro- 
magnons, an analysis of the spectral weight shows that 
long range order and antiferromagnons can account for 
only about 1/2 of the observed spectrum^]- This diffi- 
culty originates from the presence of finite spectral weight 
at high energies 0- The unidentified spectral weight in- 
dicates the presence of excitations beyond one magnon 
mode0- 

In the present paper we claim that the unidentified 
spectral weight at high energies in the spin spectrum can 
be identified with deconfined gapped spinons. In a theo- 
retical point of view, if antiferromagnetism is supposed 
to originate from algebraic spin liquid via spontaneous 
chiral symmetry breaking, the emergence of deconfined 
spinons seems to be possible. The mechanism of spinon 



deconfinement results from the existence of fermion zero 
modes in single monopole potentials. It should be noted 
that the mechanism of monopole suppression in the 
present antiferromagnet completely differs from that in 
the algebraic spin liquid, where existence of quantum 
criticality (critical fluctuations of matter fields) is the 
origin of spinon deconfinement Q. In high energy physics 
the mechanism of monopole suppression by fermion zero 
modes is well known m^SQl- I n the context of quantum 
antiferromagnetism there was a try to find fermion zero 
modes. Marston has tried to find fermion zero modes in 
the algebraic spin liquid)^. Unfortunately, the fermion 
zero mode was not found. 



In the present communication we find a fermion zero 
mode in the antiferromagnetism arising from the alge- 
braic spin liquid via spontaneous chiral symmetry break- 
ing. In the present antiferromagnet the effective action 
is pretty much similar to a well studied action in high 
energy physics, where a fermion zero mode exists. In the 
antiferromagnet Neel vectors can make a hedgehog con- 
figuration around a magnetic monopole of compact U(l) 
gauge fields. Remarkably, in the monopole-skyrmion 
composite potential the Dirac fermion is shown to have a 
zero mode. The emergence of the fermion zero mode 
forbids condensation of magnetic monopoles, resulting 
in deconfinement of Dirac spinons. Notice that Dirac 
spinons are massive owing to spontaneous chiral sym- 
metry breaking. These gapped spinons would appear to 
make broad continuum spectrum above their mass gap 
in inelastic neutron scattering^ IsL fiof . 

The main body of the present paper consists of three 
major parts. The first is to introduce a fermionic non- 
linear a model [Eq. (1)] as a low energy effective field 
theory for the quantum antiferromagnetism. The second 
is to prove the existence of fermion zero modes based on 
the proposed fermionic nonlinear a model [Eq. (3)] . The 
last is to investigate the effect of the fermion zero mode 



[Eq. (13)]. 
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II. EXISTENCE OF FERMION ZERO MODES 
AND DECONFINEMENT OF SPINONS 

A. Fermionic Nonlinear a Model for Quantum 
Antiferromagnetism 

We consider the following effective action called 
fermionic nonlinear <r model 



Z = J DnD\Dip n Da^e 
Sasl = J d 3 x 
Sm = / d 3 xm^'ip n (n ■ T nm )ip 



Sasl—Sm — <Sa 



SnLct 



M 



^IM'-'^N-'-n- c> 



In Sasl "0n is a massless Dirac spinon with a flavor index 
n = l,...,N associated with SU(N) spin symmetry, 
is a compact U(l) gauge field mediating long range inter- 
actions between Dirac spinons. e is an internal electric 
charge of the Dirac spinon. Sasl in Eq. (1) was proposed 
to be an effective field theory for one possible quantum 
paramagnetism, called algebraic spin liquid, of SU(N) 
quantum antiferromagnets on two dimensional square 
lattices^. However, the stability of the algebraic spin 
liquid has been suspected owing to instanton excitations 
of compact U(l) gauge fields [Till. Condensation of in- 
stantons (magnetic monopoles) Il2| is well known to cause 
confinement of charged particles j 131 Il4 . Il5| . here Dirac 
spinons. Recently, Hermele et al. showed that the alge- 
braic spin liquid can be stable against magnetic monopole 
excitations J3|. Ignoring the compactness of U(l) gauge 
fields a^, one can show that the Sasl has a nontrivial 
charged fixed point in two space and one time dimensions 
[(2+l)£>] in the limit of large flavors [1 [11 [l7|, identified 
with the algebraic spin liquid. Hermele et al. showed that 
the charged critical point in the case of noncompact U(l) 
gauge fields can be stable against magnetic monopole ex- 
citations of compact U(l) gauge fields in the limit of large 
flavors 3] . Condensation of monopoles can be forbidden 
at the stable charged fixed point owing to critical fluc- 
tuations of Dirac fermions. The Sasl m Eq. (1) is a 
critical field theory at the charged critical point, where 
correlation functions exhibit power law behaviors with 
anomalous critical exponents resulting from long range 
gauge interactions [lH Hs[ . This is the reason why the 
state described by the Sasl is called the algebraic spin 
liquid. In appendix A we briefly discuss how the effec- 
tive quantum electrodynamics in (2 + 1)13 (QED3), the 
Sasl in Eq. (1) can be derived from the antiferromag- 
netic Hcisenberg model on square lattices. 

However, we should remember that the criticality of 
algebraic spin liquid holds only for large flavors of criti- 
cal Dirac spinons. If the flavor number is not sufficiently 
large, the internal charge e is not screened out satisfacto- 
rily by critical Dirac fermions. Then, gauge interactions 



can make bound states of Dirac fermions, resulting in 
massive Dirac spinons. This is known to be spontaneous 
chiral symmetry breaking (S X SB)^, HJ M, El E|. 
In the case of physical SU(2) antiferromagnets it is not 
clear if the algebraic spin liquid criticality remains owing 
to the SxSB causing antiferromagnetism. It is believed 
that there exists the critical flavor number N c associ- 
ated with SxSB in the QED* [l9ll20| . But, the precise 
value of the critical number is far from consensus [2(j. If 
the critical value is larger than 2, the SxSB is expected 
to occur for the physical N = 2 case. Then, the Dirac 
fermions become massive. On the other hand, in the case 
of N c < 2 the algebraic spin liquid criticality remains 
stable against the SxSB. Experimentally, antiferromag- 
netic long range order is clearly observed in the SU(2) 
antiferromagnet. This leads us to consider the SxSB in 
the algebraic spin liquid. In Eq. (1) Sm represents the 
contribution of a fermion mass due to the SxSB. is 
a mass parameter corresponding to staggered magneti- 
zation in the context of antiferromagnetism|2lj . ft repre- 
sents fluctuations of Neel order parameter fields regarded 
as Goldstone bosons in the SxSB. f is Pauli matrix act- 
ing on the spin (flavor) space of Dirac spinons. The mass 
parameter can be determined by a self-consistent gap 
equation, given by sa e 2 exp[—2Tr/y / N c /N— 1] in the 
1/N approximation[19j. For completeness of this paper 
we briefly sketch the derivation of dynamical mass gen- 
eration in appendix B. There are additional fermion bi- 
linears connected with the Neel state by "chiral" trans- 
formations because the Sasl in Eq. (1) has more sym- 
metries than those of the Heisenberg model|22L 12^. |25| . 
These ord er p arameters are associated with valance bond 
orders[22l|2J,|25j. But, in the present paper we consider 
only the Neel order parameter in order to obtain the 0(3) 
nonlinear a model. 

Contributions of high energy spinons in Sasi^+Sm 
lead to the SnloM in the gradient expansion |23l |24| . 
This effective action is nothing but the 0(3) nonlinear 
<j model describing quantum antiferromagnetism. In the 
Snlo-m we introduced a Lagrange multiplier field A to 
impose the rigid rotor constraint \ri\ 2 = 1. In appendix 
C we give a detailed derivation of the nonlinear a model. 
The total effective action Sasl + Sm + Snl<tM called 
fermionic nonlinear a model in Eq. (1) naturally de- 
scribes the quantum antiferromagnetism resulting from 
the algebraic spin liquid via SxSB at half filling. Based 
on this fermionic 0(3) nonlinear a model we discuss 
physics of the quantum antiferromagnetism. 

In Eq. (1) we focus our attention on the Kondo- 
like spin coupling term n ■ iprip between antiferromag- 
netic spin fluctuations n and Dirac fermions ip n . This 
term shows that an antiferromagnetic excitation of spin 
1 can fractionalize into two fermionic spinons of spin 1/2. 
But, long range gauge interactions prohibit antiferromag- 
netic spin fluctuations from decaying into spinons. More 
quantitatively, massive Dirac spinons should be confined 
to form spin 1 antiferromagnetic fluctuations owing to 
the effect of magnetic monopoles of compact U(l) gauge 
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fields <z M . SxSB makes the criticality of algebraic spin 
liquid disappear, causing massive spinons. These spinons 
can generate only the Maxwell kinetic energy for the 
gauge field a M via particle-hole polarizations. It is well 
known that this Maxwell gauge theory shows confine- 
ment of charged matter fields in (2 + 1)D o wing to the 
condensation of magnetic monopoles^J 0, 1 1 5| . Spin 
1/2 massive Dirac spinons are confined to make spin 1 
antiferromagnetic fluctuations, i.e., ft = (ipfip), where 
(...) denotes a vacuum expectation value. A resulting ef- 
fective field theory for this antiferromagnet is obtained 
to be 



Zaf = / Dire 



Sty — 



d 6 x 



47T 



(IM + J _ |7?|2 J 



(2) 



where the Neel vector is given by ft = (if, n 3 ) with 
ri 3 = yT- |7?| 2 . In the last line we obtained an effec- 
tive field theory for small fluctuations of if fields around 
the Neel axis n 3 in the antiferromagnet. The if fields are 
nothing but aniferromagnons, considered to be spinon- 
antispinon composites tt^ = (ipT^ip) = nx ± in 2 with 
relativistic spectrum, u> — k in the low energy limit. The 
last term in the last line represents interactions between 
antiferromagnons. Low energy physics in this conven- 
tional quantum antiferromagnet is well described by the 
interacting antiferromagnons 26J. This is the result of 
conventional antiferromagnetism when we do not con- 
sider fermion zero modes. In the present paper we show 
that the presence of fermion zero modes can alter this ef- 
fective field theory Eq. (2) at high energies. We will see 
that the fermion zero mode can appear from the Kondo- 
like coupling term, resulting in magnon decaying into 
spinons. 



B. Existence of Fermion Zero Modes 

Now we show that a fermion zero mode can arise 
in a monopole-skyrmion composite potential. In order 
to introduce monopole potentials we separate the com- 
pact U(l) gauge field a M into a M 



a,. + a?, n , where 



a£ represent magnetic monopole (instanton) potentials 
and a^ 11 , gaussian quantum fluctuations. In addition, 
we consider skyrmion configurations ft sky . The single 
monopole and skyrmion potentials are given by a c ^ = 



a(r)e 3atJl x a and 



sky 



$>(r)xij,, respectively, where a(r) 



and <&(r) are proportional to r in r — > oo with r = 
\Jt 2 + x 2 + y 2 0, |23 ■ Integrating over the Dirac spinon 
fields in the fermionic nonlinear a model Eq. (1), we ob- 
tain the following fermion determinant in the monopole- 
skyrmion composite potential, S$ — —Nlndet 7^(3^ — 

* a /i' — i a Ji l ) + m^ft skv ■ f ■ For the time being, we ignore 



gaussian quantum fluctuations because our objective 
is to find a fermion zero mode in the single monopole 
potential a°}. In order to calculate the determinant we 
solve an equation of motion in the monopole-skyrmion 
potential 

(Tjud/Am + ia(r)(7 x x) 3 + m^^r^T^)^™ = Eip n . 

(3) 

Remember that in the absence of the skyrmion poten- 
tial (m^p = 0), i.e., in the algebraic spin liquid the 
fermion zero mode was not found|||. As mentioned 
in the introduction, the presence of the Neel vector 
makes the Dirac equation (3) have the essentially same 
structure as that utilized in the context of high energy 
physicsHHIlQ. SU(2) gauge theory in terms of mass- 
less Dirac fermions and adjoint Higgs fields interacting 
via SU(2) gauge fields has been intensively studied in the 
presence of Kondo-like isospin couplings between these 
matter fields 0, H H 0. When the Higgs fields are con- 
densed, topologically nontrivial stable excitations called 
't Hooft-Polyakov monopoles can appear |27l |28| . The 't 
Hooft-Polyakov monopoles consist of gauge-Higgs com- 
posite potentials [53, Hg. Jackiw and Rebbi showed that 
the Dirac equation in the presence of the Kondo-like 
isospin couplings has a fermion zero mode in a 't Hooft- 
Polyakov monopole potential 4] . Here, the Neel order 
parameter fields play the same role as the adjoint Higgs 
fields. 

Following Jackiw and RebbiQ, we explicitly demon- 
strate that Eq. (3) has a zero mode. We rewrite Eq. (3) 
in terms of the two component spinors xt with E = 



+iay(o- 2 )ijxf n ~ ian(?i)ijXf n ± m Tl>®%»xf m (T fiT ) 



(4) 



Inserting xf n — ■M-twTmn with a two-by-two matrix M. ± 
into the above, we obtain 

a 3 d T M ± + a 2 d x M ± + a 1 8 v M ± 

+iaya 2 M ± - iaxa\M ± T m^M ± x ll a fl = 0. (5) 

Now the spin matrices r k and the Dirac matrices a k 
are indistinguishable Finally, representing the matrix 
M ± in Mf m = g ± 8 im + g^(T% m , we obtain the coupled 
equations of motion for the numbers g ± and gj^ 

(d T =F m^T)g ± - i(d x + iay ± m^y)gf 

+i{d y — iax ± m^x)g 2 t = 0, 

(d x + iay =F m^y)g ± + i(d T ± m^$T)gf 

—i{d y — iax ± m^x)g 3 = 0, 

(d y — iax =F m^x)g ± — i(d T ± m^^T)gf 

+i(d x + iay ± m^y)gf = 0, 

(d T =F m^r)gf + (d x + iay =F m^$y)gf 

+(d y — iax =p mTf,^x)gf = 0. (6) 
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These equations yield the following zero mode equations 

(d T + m^$T)g~ = 0, 

(d x + iay + m^y)g" = 0, 

(d y — iax + m, l p<&x)g~ = 0. (7) 

A normalizable zero mode solution g~ is given by g~ ~ 



exp 



-i J dxa(r)y + i J dya(r)x 



exp 



Existence of the zero 



J dxm^(r)y — J dym^<&(r)x 

mode makes the fermion determinant zero in the single 
monopole excitation. As a result single monopole exci- 
tations are suppressed by the zero mode and thus, the 
condensation of magnetic monopoles is forbidden. The 
suppression of monopole condensation results in decon- 
finement of internally charged particles, here the Dirac 
spinons. This deconfined antiferromagnetism differs from 
the usual confined one described by Eq. (2). Antiferro- 
magnetic spin fluctuations fractionalize into spinons be- 
cause the U(l) gauge interactions are not confining any 
more. Below we discuss an effective field theory to de- 
scribe this unusual antiferromagnetism. 



C. Effect of Fermion Zero Modes: 't Hooft 
Effective Interaction and Deconfinement of Massive 
Spinons 

In high energy physics it is well known that the contri- 
bution of instantons (monopoles) 12] in the presence of a 
fermion zero mode gives rise to an effective interaction 
to fermions^, 0, |23- This interaction is usually called 
't Hooft effective interaction. In order to obtain the ef- 
fective fermion interaction it is necessary to average the 
partition function in Eq. (1) over various instanton and 
anti-instanton configurations. Following Ref. j2^|, we 
first consider a partition function in a single instanton 
potential 

V'&n} = f d 3 x(Mx) ltl dMxf 



(8) 



Here $^ is the zero mode obtained from Eq. (7). A 
fermion mass m is introduced. Later the chiral limit 
m — ► will be taken. The effective action including 
the effective potential [?/>„] in Eq. (8) gives a cor- 
rect green function in a single instanton potential|23|. 

S'(z,y) =< M*)My) >= - * i(x) * Uv) +S (x,y) with 
the bare propagator So(x,y) = ('y fJ _d ll )~ 1 d(x — y). Thus 
the partition function in Eq. (8) can be used for average 
of instantons |23j. The partition function in the presence 
of N+ instantons and iV_ anti-instantons can be easily 



built up [21 



Difj n Da q u u e- 5 iS ^ft-0*» 



(m- < V 1 ^] >Y + (m- < V 1 ^} >Y . (9) 



qu 



Here we admit the noncompact U(l) gauge field 
representing gaussian quantum fluctuations. In the fol- 
lowing the index qu is omitted. (...) means averaging 
over individual instantons. Introducing instanton aver- 
aged nonlocal fermion vertices Y± = —V < V [ip n ] >— 
— J dPzjrftV 1 ^' [ipn] with volume V, where Zj/j\ repre- 
sent instanton (anti-instanton) positions [2^, we obtain 
the following partition function in the chiral limit m — > 

Z^= I D-iP n Da tl e- f dax ^~<^- ia ^ 



d\ ± 



dT±e 



iA+(y+-r + )+Jv + in ^.+(+-- ) 



(10) 



Integration over X± and T± recovers Eq. (9) in the chi- 
ral limit. In the thermodynamic limit N±, V — > oo and 
N±/V fixed, integration over r± and A± can be per- 
formed by the saddle point methodic . Integrating over 
r± first, we obtain 



dX± A f +(ln^-l)+(+^-) 

2tt 6 



Dip n Da^e~ I d3x ^ ( - d »- ia ^ +iX + Y++iX - Y -(ll) 



An explicit calculation for the instanton av erag e shows 
that the vertex y ± corresponds to a mass[23j, = 



Here 



/ (f^ (2TrpF(k)) 2 {p n ^^ n with 75 = , 7 () 

F(k) is associated with the fermion zero mode in the 
effective potential y I [V' n ] in Eq. (8). In the present 
paper we do not perform an explicit calculation for the 
instanton average in Y± and thus, we do not know the 
exact form of F(k). Our objective is to see how the 't 
Hooft interaction appears as an instanton effect. Here 
p is the size of an instanton. Although the instanton 
(magnetic monopole) can be considered to be a point 
particle, the instanton-skyrmion composite would have 
its characteristic size p. In order to determine the size 
p we should solve an equation of motion for the Neel 
vector ft in the presence of an instanton configuration. 
In the present paper we do not examine this issue. We 
assume its existence. Owing to the neutrality condition 
of magnetic charges, N + = N- = N]/2 is obtained in 
Eq. (11), where Nj is the total number of instantons 
and anti-instantons. The saddle point solution of A + = 
A_ = A in Eq. (11) gives rise to cancellation of the 75 
term in the mass, causing a momentum dependent mass 
m(k) = m I F 2 {k) with m 7 = A(2?rp) 2 23]. The mass mj 
is determined by the saddle point equation for A usually 
called a self-consistent gap eauation|23| 



d 3 k 



\k) 



N T /V J (2tt) 3 k 2 +m 2 (k) 



(12) 
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Ignoring the momentum dependence by setting F(k) = 
1 for simplicity, we obtain the 't Hooft mass mj — 
( \ 1/2 

2A T/2 ( ^ ) with a momentum cut-off A in the small 
mass limit. Since the mean density of instantons is 
proportional to the instanton fugacity, Ni/V ~ Vm = 
e -S in3t w jth an instanton action Si nst ~ l/e 2 [l3l fl5| . 

1/2 

the fermion mass is roughly given by mi ~ ym ■ We 
can obtain the following effective Lagrangian in terms 
of Dirac spinons ip n with the 't Hooft effective mass 
mi interacting via noncompact U(l) gauge fields a M , 
£</j = tpnJuidf, - idfjipn + mi^ntpn- The formal appear- 
ance of the fermion mass does not necessarily mean that 
the fermions are massive. This is because the fermion 
mass is determined by density of instantons. In order to 
determine the fermion mass, we should find the instanton 
fugacity, i.e., reveal the state of instantons. 

It is one of the most difficult problems to determine 
the state of monopoles (instantons) in the presence of 
matter fields. In this paper we do not pursue to deter- 
mine the state of monopoles precisely. Instead, we con- 
sider all possibilities. Generally speaking, there would 
be three possible monopole states. These are monopole 
plasma, monopole dipolar and monopole "liquid" phases. 
These phases are characterized by y m — > oo, y m — > and 
y m 7^ 0, respectively, where y m is the monopole fugacity. 
We exclude the first, monopole plasma phase because the 
presence of fermion zero modes does not allow the con- 
densation of monopoles. 

At first glance a liquid phase of magnetic monopoles 
sounds quite strange. It is known that the liquid state 
does not a ppea r in the Abelian Higgs model without 
fermions [29t f30L l3ll l3^| . But, in the present case we do 
not have any evidence to exclude this monopole state._ In 
(2 + 1)D the basic trend is confinement, i.e. , 
Il5| away from quantum criticalityji 
ing to the confinement tendency many instantons would 
be excited although the presence of fermion zero modes 
does not allow instanton condensation. Thus, we should 
treat dense uncondensed monopoles. We claim that in 
order to solve this problem a new methodology beyond 
the dilute approximation of monopoles is required. As 
far as we know, this methodology is not found yet. Note 
that the usual renormalization group equation for the 
monopole fugacity is obtained from the dilute approxi- 
mation of monopolesH ES [H |H El IH . It is not 
clear if this standard renormalization group equation is 
applicable to high density limit. In this respect we have 
no clear evidence to exclude the monopole liquid phase 
in the presence of fermion zero modes. We may view the 
emergence of the liquid state as the proximate effect of 
the Higgs-confinement phase in the presence of fermion 
zero modes. There exist some reports about a new phase 
instead of plasma and dipolar phases in two dimensional 
Coulomb gas when the density of particles is high[34j. 
Furthermore, a new fixed point with nonzero monopole 
fugacity was recently reported even in the QED3 only 
with massless Dirac fermions [l7j. 



mis monopoie state, m 
merit, i.e., y™ — ► 00 [TH 

in m mm. ow- 



In the liquid phase of monopoles Dirac spinons obtain 
the 't Hooft effective mass (mi) while in the dipolar phase 
of monopoles the 't Hooft mass vanishes. But, it should 
be noted that we are considering antiferromagnetism. In 
antifcrromagnetism Dirac spinons are massive (m^) via 
SxSB. This mass completely differs from the 't Hooft 
mass. Although it is not easy to evaluate the 't Hooft 
mass, we expect that the chiral mass m^ would be larger 
than the 't Hooft mass mj. This assumption is based 
on the fact that the chiral mass m^ would be identified 
with the staggered magnetization observed in quantum 
antifcrromagnets. In both phases of magnetic monopoles 
the deconfined Dirac spinons would be massive owing to 
the antiferromagnetism. At present we do not have any 
idea how to distinguish these two monopole states. 

Combining the existence of fermion zero modes in sec- 
tion II-B and their physical effects in section II-C, we can 
reach the following effective Lagrangian from Eq. (1) in 
the deconfined antiferromagnet 

Zaf = J D^Da^Dne- 1 d3xLAP , 
£af = ipnlfiidfi - ia^ipn + m^ n Tl m il) m + —z\d x a\ 2 



^3 „/. \~\2 



-^m^ij) n T nm ij) m \ii\ + m^ip n T nm ip m ■ tt 



-in 



(\d^\ 2 + (n-d^) 2 ), 



(13) 



where we used n 3 = y/l - \tt\ 2 « 1 - \\tt\ 2 as Eq. (2). 
In this effective Lagrangian we should understand that 
the gauge field a M is noncompact in Eq. (13) owing to 
the effect of fermion zero modes while the gauge field 
in Eq. (1) is compact, causing confinement of Dirac 
spinons and thus, resulting in Eq. (2). In this anti- 
ferromagnet low lying excitations are antiferromagnetic 
spin waves and (noncompact) U(l) gauge fluctuations. 
Both gapless excitations cause C v ~ T 2 , where C v is 
specific heat and T, temperature 21] . Furthermore, the 
gapped Dirac spinons would be deconfined to emerge 
above their mass gap. These gapped spinons interact 
with the antiferromagnons as shown by the fourth and 
fifth terms (3^. In Fig. 1 we show schematic spin spec- 
trum in inelastic neutron scattering. Here x" (lu) is the 
imaginary part of the transverse dynamic spin suscepti- 
bility at momentum (tt,tt), given by x"( tL ') ~ X"tt(w) + 
+ x"c(w) where x\ ~ (n + ir~) is a magnon sus- 
ceptibility, x"^ ~ (iJjt + ip^r~ ip) , a spinon susceptibility, 
and x"c ~ (■0r + '07r _ ) + (tt + %\)t~ if>) , a "coupling" sus- 
ceptibility in a highly schematic form. At low energies 
antiferromagnons give dominant spectral weight (x" v) 
while at high energies Dirac spinons exhibit broad con- 
tinuum spectrum (x" i>)- In our scenario the unexplained 
1/2 spectral weight in the inelastic neutron scattering^ 
would be identified with the gapped Dirac spinons. Al- 
though scattering between spinons and magnons in Eq. 
(13) does not alter the spin spectrum in Fig. 1 quali- 
tatively, elaborate calculations are required in order to 
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FIG. 1: Schematic spin spectrum in inelastic neutron scat- 
tering 



understand the spin spectrum more quantitatively and 
precisely. 



III. SUMMARY 

In summary, we investigated the two dimensional 
quantum antiferromagnet resulting from the algebraic 
spin liquid via spontaneous chiral symmetry breaking. 
This antiferromagnet is described by the fermionic non- 
linear a model, Eq. (1) in terms of Dirac spinons in- 
teracting via not only U(l) gauge fluctuations but also 
antifcrromagnetic spin fluctuations n. We showed that 
the Kondo-like spin couplings between the Dirac spinons 
and Neel vectors give rise to the fermion zero mode in the 
single monopole potential. The existence of fermion zero 
modes suppresses the condensation of monopoles, thus 
causing the deconfinement of spinons. As a result we ob- 
tained the effective field theory Eq. (13) in the deconfined 
antiferromagnet, which differs from the conventional con- 
fined antiferromagnet described by the 0(3) nonlinear a 
model Eq. (2). From the effective field theory Eq. (13) 
we argued that the deconfined massive spinons would be 
observed as particle-hole continuum above their mass gap 
in the dynamic spin susceptibility [3. 1^. ITo) . 
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APPENDIX A 

In appendix A we briefly sketch how we can ob- 
tain the Sasl in Eq. (1) from the antiferromagnetic 
Heisenberg model on two dimensional square lattices, 



<»,i> 



Si ■ Sj with J > 0. Inserting the follow- 



ing spinon representation for spin, Si = \f\ a T a jifiji into 
the above Heisenberg model, and performing the stan- 
dard Hubbard-Stratonovich transformation for an ex- 
change interaction channel, we obtain an effective one 
body Hamiltonian for fermions coupled to an order pa- 
rameter, H ef f = -J E<i, 3 > flaXijfja ~ h.c. Here f ia 
is a fermionic spinon with spin a =t, |, and x.ij is an 
auxiliary field called a hopping order parameter. No- 
tice that the hopping order parameter xij is a com- 
plex number defined on links if. Thus, it can be de- 



composed into Xij 



I Xij I 



where 



Xij 



and 0a are 



lowing U(l 

0'. - a ' 

ij 



the amplitude and phase of the hopping order param- 
eter, respectively. Inserting this representation for the 
Xij into the effective Hamiltonian, we obtain H e ff = 

_ ^Z)<i,i> \Xij\fia e%e<i fia ~ h.c. We can easily see that 
this effective Hamiltonian has an internal U(l) gauge 
symmetry, H' eff [f' ia , 6'^] = ff e //[/ia,%] under the fol- 
phase transformations, f' ia — e l<l>i f ia and 
— cf>j . This implies that the phase field fty 
of the hopping order parameter plays the same role as a 
U(l) gauge field ay. When a spinon hops on lattices, it 
obtains an Aharnov-Bohm phase owing to the U(l) gauge 
field dij . It is known that a stable mean field phase is a tt 
flux state at half filling [Hl2l|. This means that a spinon 
gains the phase of it when it turns around one plaquette. 
In the 7r flux phase low energy elementary excitations are 
massless Dirac spinons near nodal points showing gap- 
less Dirac spectrum and U(l) gauge fluctuations [H |21|. 
In the low energy limit the amplitude \xij\ is frozen to 

\Xij\ = -^1 < fjafia > |. A resulting effective field the- 
ory for one possible quantum disordered paramagnetism 
of the antiferromagnetic Heisenberg model is QED3 in 
terms of massless Dirac spinons interacting via compact 
U(l) gauge fields, Sasl in Eq. (1). 



In the Sasl tp n 



_ ( xt 

Xn 



is a four component mass- 



less Dirac fermion, where n = 1,2 represents its SU(2) 
spin (t, I) and ± denote the nodal points of (7r/2, ±7r/2) 
in momentum space. Usually, SU(N) quantum antiferro- 
magnets are considered by generalizing the spin compo- 
nents n = 1,2 into n — 1,2, ...,N. The two component 



spinors Xn are given by xt 
and X2 = 



X 2 = 



fl2o 



/T2o 
/t2e 



Xl = 

, respectively. In 



file _ ^ _ 

/llo J ' " A2 V fi 2e 

the spinon field f a bc a =T,I represents its SU(2) spin, 
b = 1,2, the nodal points (+,—), and c = e, o, even 
and odd sites, respectively [2l|. Dirac matrices 7^ are 



given by 70 = 



0-3 









/ 0-2 





-0-3 / 


) 71 = 


{ 


-0-2 J 



and 



72 



o~ 1 

-(71 

Clifford algebra [7^,7^] + 



, respectively, where they satisfy the 
2V gll. 
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APPENDIX B 

In appendix B, for completeness of this paper we briefly 
sketch how we obtain the dynamically generated spinon 
mass in Sm in Eq. (1). A single spinon prop- 
agator is given by G~ 1 (k) — G " 1 (fc) — £(fc), where 
@o s ) = *7m^m ^ s * ne hiverse of a bare spinon prop- 
agator, and £(&;), a spinon self-energy resulting from 
long range gauge interactions. The spinon self-energy is 
determined by the self-consistent gap equation, £(£:) = 

J j^ Tr h^ G ( k -<lhv D ^( < l)}i where D^(q) is a renor- 
malized propagator of the U(l) gauge field a M due to 
particle-hole excitations of massless Dirac fermions. The 
self-energy can be written as = — m^(k)r 3 for stag- 
gered magnetizationj^. Inserting this representation 
into the above self-consistent gap equation, we obtain 
the following expression for the spinon mass, m^(p) = 

I fcCm^ ^VCP - k )- Thc renormalized gauge 



The second term in the last line is well derived in Ref. 



— In det 
2 



1 



-d 2 + m?,, 



N, 



Tr / d 3 x(x\ln 



d 3 k 



-d 2 



c Trln 



4 ' ^ 



(2tt) 3 



-<9 2 + rr& 
d 3 k 



Trh 



(2tt) 3 
m^^d^n ■ f) 



k 2 + m 2 - 2*^ - 9 2 



propagator D )lu (q) is obtained to be D flll (q) w II J (q) 
f^jui/— ^r 1 ) in the Lorentz gauge [2l|. where n /iI/ (g) 



A 



d 3 fc (m^^d^n ■ f) 



— Tr( 

Tr) 3 V 



47T 



(27T) 



fc 2 



(C2) 



-A / ^p-Tr[7 Al Go(fc)7^Go(fc - q)] is the polarization 

function of massless Dirac fermions in the 1/A approx- 
imation. Inserting this gauge propagator into the gap 
equation and performing an angular integration, one can 



In the above Tr stands for not a functional but a usual 
matrix trace for both flavor (spin) and spinor indices. In 
going from the third to the fourth line we have dragged 
the factor e lkx through the operator, thus shifting all dif- 



find (p) = 



4 



Jo dk 



'■' r '^( fc ) (]„ 



■p— \k — p\), where ferential operators d : 



ik^ |2_ 



Expanding the 



A is a momentum cutoff 21]. This integral expression is 

2 dm^jp) \ _ 



dp 



equivalent to the differential equation, (p 
_ 8 J^4£L with boundary conditions, K^btM. + 

w'N p 2 +m^(p) J ' dp p= A 

m^(A) = and < m^(0) < oo[H|2l|. In Ref ' G3 
this equation is well analyzed in detail. Its solution is 
given by my w e 2 exp[—2ir/y/N c /N — 1] in the case of 
N<N c \i]l\2^. 



APPENDIX C 

In appendix C we discuss how the 0(3) nonlinear a 
model can be derived from the effective action Sasl + 
Sm in Eq. (1). Integration over the Dirac fermions in 
Sasl + Sm results in the following expression 



S 



NLcrM 



= -lr 



Dtp n exp{- 



d 3 :, 



' 7~ nrn 



argument of the logarithmic term in powers of d^ft and 
of 2ik fJ ,d fl + d 2 , one can easily obtain the expression in 
the fifth line. Performing the momentum integration, 
we obtain an effective spin stiffness proportional to the 
mass parameter m^. This implies that the rigidity of 
fluctuations in the Neel field is controlled by the mass 
parameter of Dirac spinons. Eq. (C2) is nothing but 
the 0(3) nonlinear a model describing quantum antifer- 
romagnetism. Note that higher order derivative terms in 
the gradient expansion are irrelevant in (2 + l)D in the 
renormalization group sense. 

Next, we sketch the derivation of the Maxwell gauge 
action. Expanding the argument of the logarithmic 
term in Eq. (CI) to the second order of the gauge 

field a M , we obtain S gauge = j -^ s \a ll {q)Jl lxv {q)a v {-q), 
where the fermion polarization function IL. t „ (q) is 
given by U^q) = -N J ^Tr[^G(k + qh»G(k)} 
with the single spinon propagator G(k) = [ij^k^ + 
m^n ■ f\ . Utilizing the Feynman identity and trace 
identity of Dirac gamma matrices, one can obtain 
Unl {d ~ ia )i> n ^ ne f°ll° w ing expression for the polarization function, 
- ' " '' " H uv (q) 



x)x(m 2 h +q 2 x(l 



-N In det 
A 



2 

d 3 x 



2 i 2 

Am^,„ Jl0 A 



In det{-& 2 + m^) 



-4m:. 



2N(TrI)^^(q% v - q»q v ) jj dx(l 
This 



■'( 



211 



leads to the 



4tt 



l^n] 



\d x a\ 



(CI) 



4>J 3 V a/4 

Maxwell gauge action in Eq. (CI). 

Lastly, we should comment the reason why imaginary 
terms do not arise in the present nonlinear a model be- 
cause some previous studies have shown the emergence 



8 



of imagi nary terms[36|, |3?J. Following the evaluations 
in Ref. |36|. one can obtain two imaginary terms in 
the irreducible representation of gamma matrices; one 
is a coupling term ia M J M between topologically nontriv- 
ial fermionic currents J M = ■^e^ v xe a ^ 1 n a d v n^ dxn 1 and 
U(l) gauge fields a M , and the other, a geometrical phase 
^-/V7rr[n]|36l l37j. The key point is the representation of 
Dirac gamma matrices\3a(. Here, we utilized four-by- 
four Dirac matrices by combining the two nodal points. 
Note that the signs of Pauli matrices in the Dirac gamma 
matrices are opposite for the nodal points ±. This fact 
results in cancellation of the imaginary terms. The first 
imaginary term can be considered from a variation of the 
logarithmic term in Eq. (CI) with respect to the gauge 
field a,. 



is also exactly zero owing to the same reason. Another 
way to say this is that the signs of mass terms for the 
Dirac fermions (x^ an d Xn) a ^ ^e two Dirac nodes (+ 
and —) are opposite, resulting in cancellation of the parity 
anomalu\3?fy . If we fix the Neel vector in the z direction 
(n = z), we can see the opposite signs explicitly from 

m ipXn' J z' r n mXm = m "4>Xn T nmXm 

massive Dirac fermions (xt) contribute to the imaginary 
terms, respectively. However, the signs of the imaginary 
terms are opposite and thus, the cancellation occurs. As 
a result the imaginary terms do not appear in the present 
nonlinear a model. This was already discussed in Ref. 



~- rn^xta z T^ m xtr - 

m^Xn^rnXm- B Oth 



iNTr 
= iNTr 
iNTr 



1 



M 7m - ia^) + ™</> (n-T)- 



8 a 



ia^)] 2 



>^ + m J — 



-d 2 



(C3) 



In this expression the key point is the triple product of 
Dirac matrices, 7^7j,7a in the last line. In the irreducible 
representation of gamma matrices, i.e., Pauli matrices, 
this contribution is nonzero, leading to t nv x- As a result 
the ima gina ry term of ia ll e lllJ \£ a p~ f n a d u nP dxn 1 can be 
obtained |36ll37| . Another e Q( 3 7 associated with the Neel 
vectors appears from the triple product of Pauli matrices, 
T a TpTy. On the other hand, in the present representation 
of Dirac matrices the contribution of the + nodal point 
leads to +e nu x while that of the — nodal point, —e^ v x- 
Thus, these two contributions are exactly cancelled. The 
geometrical phase term, considered from a variation of 
the logarithmic term in Eq. (CI) with respect to the 
Neel field nlMEl, 



-NTr 



-NlmTr 



liidf, + m i ,(n-T). 



NTr 



m^8{n ■ t)- 



-d 2 



-d 2 + m 2 , - m^^d^n ■ f) 
- ml V -d 2 + m 2 J r 



It is possible that the mass terms have the same 
signs. Considering the two gamma matrices of 74 

1\ , for 

and 75 = 



-I 



I21l. we can obtain the 



terms with the 



same 

x + - 



nm A-m 



signs, C M 

m ipXn T nmXr, 



I 

following mass 

in 

These mass terms can arise from the algebraic spin liq- 
uid, Sasl m Eq. (1) via SxSB because the algebraic 
spin liquid has the enlarged symmetry [3 S3- as dis- 
cussed earlier. In this case the cancellation dose not oc- 
cur and thus, the imaginary terms necessarily arise. This 
antiferromagnetism would not be conventional since it 
breaks not only time reversal symmetry but also parity 
symmetry. When this antiferromagnetism disappears via 
strong quantum fluctuations, its corresponding quantum 
disordered paramagnet is expected to be the chiral spin 
liauid\24. 139L |40| . In the present paper we did not discuss 
the chiral spin liquid. 
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